High-temperature surface superconductivity in topological flat-band systems 
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We show that the topologically protected flat band emerging on a surface of a nodal fermionic 
system promotes the surface superconductivity due to an infinitely large density of states associated 
with the flat band. The critical temperature depends linearly on the pairing interaction and can 
be thus considerably higher than the exponentially small bulk critical temperature. We discuss an 
example of surface superconductivity in multilayered graphene with rhombohedral stacking. 
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Normal Fermi liquid is the generic form of a system of 
interacting fermions. Fermi liquid has a finite density of 
states (DOS) at zero energy, which may lead to instabil- 
ities at low T with formation of broken symmetry states 
with smaller DOS. However, there is a class of fermionic 
systems with diverging DOS: the systems with a disper- 
sionless spectrum that has exactly zero energy, i.e., the 
flat band. Historically this was first discussed in the con- 
text of Landau levels. However, flat bands may emerge 
also without a magnetic field, for example in strongly 
interacting condensed matter systems^—, in layered sys- 
tems with integer-valued pseudospin~, in 2+1 dimen- 
sional quantum field theory dual to a gravitational the- 
ory in the anti-de Sitter background 6 , etc. In some cases 
the flat band is protected by topology in the momen- 
tum space: topologically protected zero modes emerge 
in cores of quantized vortices^—, on surfaces of gapless 
topological media such as nodal superconductors 9 - - — and 
multilayered graphen o 9 ' 12 " — , as well as at the edges of 
graphene sheets 9 -^. 

In this report we consider a three dimensional (3D) sys- 
tem where the topologically protected flat band with its 
singular DOS appears on the surface giving rise to the 2D 
surface superconductivity. This property is generic and 
does not depend much on the details of the system. For 
illustration we use the multilayered graphene with rhom- 
bohedral stacking, where a surface flat band appears in 
the limit of large number of layers. We show that the su- 
perconducting critical temperature depends linearly on 
the pairing interaction strength and can be thus consid- 
erably higher than the usual exponentially small critical 
temperature in the bulk. This may open a new route to 
room-temperature superconductivity 9 . Formation of sur- 
face superconductivity is enhanced already for a system 
having N > 3 layers where the normal-state spectrum 
has a power-law dispersion £ p cx \p\ N as a function of 
the in-plane momentum p. The DOS v(£ p ) oc £p 2 N ^ N 
has a singularity at zero energy which results in a dras- 
tic enhancement of the critical temperature. We also 
demonstrate that doping leads to a suppression of the 
surface critical temperature, contrary to its effect on the 
bulk superconductivity where the critical temperature is 
increased— li 6 -. 



a. The model. We consider multilayered graphene 
structure of N layers in the discrete representation with 
respect to the interlayer coupling. We choose the rhom- 
bohedral stacking configuration considered in&i^r— and 
assume for simplicity that the most important are jumps 
between the atoms belonging to different sublattices pa- 
rameterized by a single hopping energy t. More gen- 
eral form of the multilayered Hamiltonian can be found 
in Refsi 17 ' 18 . In the superconducting case the Hamilto- 
nian has the form of a matrix in the Nambu space. The 
Bogoliubov-de Gennes (BdG) equations are 
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where the sum runs over the layers. The normal-state 
Hamiltonian^ 

Hj = v F (a ■ p)S l . J - ta+Sij+i - taS^j-i , (1) 

& = {&x-> %)j &± — (^ x zLia y )/2, and Vi are matrices 
and spinors in the pseudo-spin space associated with two 
sublattices. This Hamiltonian acts on the envelope func- 
tion of the in-plane momentum p taken near one of the 
Dirac points, i.e., for |p| <C h/a where a is the interatomic 
distance within a layer; vf = 3toa/2h where to is the the 
hopping energy between nearest-neighbor atoms belong- 
ing to different sublattices on a layer. The particle-like, 
iii, and hole- like, wave functions near the Dirac point 
are coupled via the superconducting order parameter Aj 
that can appear in the presence of a pairing interaction. 
Here we do not specify the nature of pairing which can be 
either due to the electron-phonon interaction or due to 
other interactions that have been suggested as a source 
for intrinsic superconductivity in graphene, see Refs. 19 . 
As a reasonable starting point we assume s-wave sym- 
metry of the order parameter and neglect fluctuations 
for simplicity, though they could, in principle, be rele- 
vant for 2D superconductivity. The excitation energy for 
particles and holes is measured upwards or downwards, 
respectively, from the Fermi level which can be shifted 
with respect to the Dirac point. We assume that the 
shifts at the outermost layers may be different from the 
bulk chemical potential due to the presence of a surface 
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charge, i.e., fa = fi for i ^ 0, N while /ii,jv = A* + &Hi,n- 
The order parameter and the Fermi level shifts fa are 
scalars in the pseudo-spin space. We assume that A; 
and fa are much smaller than the inter-layer coupling 
energy t > 0, which in turn is t <C £o- Usually, t ~ 0.1 to 
where to ~ 3 eV— . 

b. Spectrum. We decompose the wave function 



or 
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into the spinor functions localized at each sublattice 
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The BdG equations take the form 



t 3 [vf(Px - ip y )a n - - p,&+] = Ea+, n^l, (3) 

t 3 [u^Ox + ^)a£ - ta+ +1 - na~] = E&~, n + N. (4) 
We introduce matrices and vectors in the Nambu space 
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In Eqs. © and Q we assume that A„ 7^ only at the 
outermost layers, while A„ = for n ^ 1,7V. The argu- 
ments supporting this model are given below. We also 
neglect A„ as compared to t in Eqs. ^ and ((H) for n = N 
and n = 1, respectively, as they lead to higher-order cor- 
rections in A/i. The particle and hole channels are thus 
decoupled if n 7^ 1,N which determines the coefficients 
= A ± e lp ^ dn and the energy in terms of the transverse 
momentum p z (d is the interlayer distance}^ 
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where p = ^Jpl+p^ and e 1 ^ = (p x + ip y )/p. 

A finite order parameter A couples the particle and 
hole channels at the outermost layers, i — 1 and i = N, 



f 3 v F (p x - ip y )a 1 -f 3 fadf = Eaf - Axaf , (6) 
f 3 v F (p x + ip y )a+ - f 3 faya N = Ea N - A N a~ . (7) 

Boundary conditions ([6]), (J7J) select p z and determine 
27V particle and hole branches of the energy spectrum. 
Looking for the branches that belong to the surface states 
with energies of the order of A and fa we solve these 
equations for E -C t. Since Eqs. (|3]), dU do not contain 
A, one can use the coefficients as obtained in Ref.— 



C_ 
71 



(irY A " 



/V F p\X-n v F p(f 3 E + (J,) 



\ t J vpp 2 — t 2 



A~ 



C 
V2 



^v F p^ N ~ n 



r vpp y 1 -- 1 v F p(f 3 E + fj) 



V t ) 



lv 2 - 1 2 



A+ 



vpp 



i{n-\- 



a i(n-f )0 



Here C is a normalization constant. We include the first- 
order corrections in energy. Having an imaginary mo- 
mentum p z for v F p < t, these solutions decay away from 
the surfaces and thus they describe the surface states. 
The vectors A ± = (A ± , B^) 71 do not depend on n. 
Equations ([6]) and (jTj) yield 
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f 3 £ p A- = (E + f 3 fa)A+ - A X A+ , 
f 3 £ P A+ = (E + f 3 fi N )A- - A N A~ , 

where £ p = t (v F p/t) , fa,N = A + df-i.N, and 

E = E(1- vpp 2 jt 2 ), // = A(l - v%p 2 /t 2 ) . (10) 
Equations 1(5)), © provide the surface-state spectrum 
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If Ai = An we have from Eq. (fTTj) 
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where S/J, — (fa — fav)/2 and jig = (fa + jj,jsr)/2. Equa- 
tions ©, © determine four independent states. If 
/i = Sfx = they are (i) E\ — E and Af = u, Bf = v, 



(ii) E 2 



-E and Ai 



v, B? 



-u, (iii) E 3 = E 



and Af = ±v, Bf — ±u, (iv) E± = —E and Af = ±u, 
Bf = qpt>. Here E = 
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The overall normalization requires ^X^n^iO^n I 2 
\P~ | 2 ] = 1. For £ p < t this gives 
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FIG. 1: (Color online): Spectrum of surface states for differ- 
ent numbers of layers. Left: N — 5 and right: iV = 50. The 
symbols have been calculated by exact diagonalization and 
solid lines are computed from Eq. (|12p up to the point where 
the approximation in it is valid. The three cases are: nor- 
mal case with fi n = (blue circles), A = 0.05t, /i n = (red 
squares), and A = 0.05t, fa = /ijv = 0.03t (magenta crosses). 



3 



Note that Eqs. ([TT )l -([i"4" ]) hold for £ p < t. The spectrum 
is plotted in Fig. [JJ 

If N — > oo and £ p —> for any VFp/t < 1, the surface- 
state part localized at n — N (with the coefficients A~) 
decouples from that (with A + ) which is localized at n = 
1. For a "flat band" f p 0, Eq. ([Til) y ields 



E\ = fi* N + \A N \ 2 or E 2 _ = ft + | Ai| 



(15) 



This shows that the definite signs in Eq. (fT2|) belong to 
the surface states localized at the corresponding layers, 
c. Flat band; zero doping. The gap at a layer N is 



x[l-2/(£ p , fc )], 

where f(E) is the Fermi distribution function. We as- 
sume that the cut-off momentum p c of the pairing po- 
tential V is larger than pfb = tjvp. The sum includes 
one n = N surface state which we label by k = and the 
bulk states specified by the transverse momenta p z (k), 
where k = 1, 2, . . . , N — 1 with the spectrum of Eq. ([5]). 
Therefore, Ajy = As + As where the surface contribu- 
tion comes from the flat band area p < pfb , 
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The bulk contribution comes from the momenta p > pfb • 
For such momenta, the surface state k — will also ex- 
tend to the bulk giving rise to (for T = 0) 
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+ ^Tr[w A r(p,fc)^(p,fc)] . (17) 
fc=i J 



All the bulk states with p > pfb are normalized to the 
sample width W = dN, i.e., u*(z) ~ 1/y/W. According 
to Eq. ©, E - v F p > t in Eq. fl7). Therefore, 
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If there was only the bulk contribution (A = A# = Ajy), 
the gap equation would have a nonzero solution only for 
a potential strength higher than a certain critical value 
Vp c /4irh 2 VFd > 1, as is the case in the usual single-layer 
graphene with zero dopin g 15 ' 16 . 

The surface states for p < pfb are normalized accord- 
ing to Eq. ([U]). We find from Eq. ([TB]) 
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For simplicity we assume that V is constant up to the 
cut-off momentum p c . Here u and v are determined by 
Eqs. (HH), (HH). In the case of a fiat band uv — 1/2 while 
E = A(l - v 2 F p 2 /t 2 ). For T = it gives 



A s = A = . 9 /8vr 



(19) 



where g = V^g/ft 2 is the characteristic pairing energy, 
V = V/d is the two-dimensional pairing potential. 

The ratio of the order parameter in the bulk to that 
on surface is of the order (A/t)(vFPc/t). Since A -C t, 
the contribution from the bulk states with E > t can be 
neglected if the cut-off momentum of the interaction p c 
does not considerably exceed t/vp- We thus arrive at the 
central result of our paper, namely that the surface su- 
perconductivity in the presence of a flat band dominates 
over the bulk superconductivity. This follows from an in- 
finitely large density of states associated with the flat 
band. The critical temperature is determined by Eq. (fT9)) 
with A — > 0, which gives A = 3fcsT c . Due to its linear 
dependence on the interaction strength, the critical tem- 
perature is proportional to the area of the flat band and 
can be essentially higher than that in the bulk. 

For a flat band £ p = with p c = pfb the only char- 
acteristic values in the superconducting surface state are 
the energy A and the momentum pfb- Therefore, the 
coherence length should be of the order of the only avail- 
able length scale, £o ~ ^/.Pfb- It is much larger than the 
interatomic distance, £o 3* a i since pfb <C Po ~ h/a. 

Doping destroys the surface superconductivity. This 
can be seen from Eq. (TT5|) with uv = A/2E + and E = 
(1 - v 2 F p 2 /t 2 )E+ where E+ is taken from Eq. ([UJ. The 
critical temperature is found by putting A = 0. For 
example, if \i and /kjv have the same sign, both Ao and 
T c vanish at the critical doping level that satisfies 
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If [in = n the critical doping is =2ksEc- 

d. Surface superconductivity in a finite array. Since 
the normal-state DOS defined as 
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has a low-energy singularity for N > 2, the surface su- 
perconductivity is favorable already for a system with 
a finite number of layers N > 3. A simple expres- 
sion for the zero-temperature gap can be obtained if 
N > 5. For a finite N, the value £ p can reach values 
larger than A. We use Eqs. (fT2j) - (fl4l) for zero doping 
in Eq. (fT8|) where the upper limit of integration p c is now 
such that £ c = t(vFPc/t) N ^ A. Transforming to the 
energy integral with the normal-state DOS Eq. (I2TH) we 
see that, for N > 4, the integral converges at £ p ~ A or 
P ~ Pa — PFB{A/t)7r . The zero-temperature gap is 



A = t (JL 
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FIG. 2: (Color online) Left panel: Zero-temperature gap as 
a function of the momentum cutoff p c for various N (solid 
lines). The gap saturates at p c ~ pa and approaches Eq. (|19[) 
for iV — » co. The dashed lines show the dispersion £ p for each 
N. Right panel: the self-consistently calculated A(z) profile 
at different layers, z — nd. On both panels g = O.Olt. 
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For JV > 1 we have oln = 1. The flat-band re- 
sult, Eq. (|19p. is recovered if the number of layers is 
N 3> 2 ln(i/A ). The coherence length for a finite system 
is £o ~ ti/pA- It approaches H/pfb for TV — » oo. 

The gap obtained by numerical integration of Eq. (|18[) 
with a cut-off p c is plotted in Fig. [2] left panel. The right 
panel of Fig. [2] shows the order parameter as a function 
of the transverse coordinate. It extends into the bulk 
only over a few interlayer distances due to a decay of the 
wave functions. Taking this into account we have chosen 
the model, Eqs. (|3J— in which the order parameter is 



nonzero only on the outermost layers. 

e. Conclusion. The flat band with infinite DOS 
emerges in semi-metals with topologically protected 
nodal lines. The flat band promotes surface supercon- 
ductivity with T c proportional to the pairing interaction 
strength and to the area of the flat band in the momen- 
tum space which is determined by the projection of the 
nodal line onto the surface. The critical temperature can 
thus be considerably higher than the exponentially small 
T c in the bulk. Formation of surface superconductivity 
is enhanced already for a system with a number of layers 
N > 3 where the normal DOS has a singularity at zero 
energy. Topologically protected flat bands may also ap- 
pear on interfaces, twin boundaries and grain boundaries 
in bulk 3D topological materials leading to an enhanced 
bulk T c . Indications towards surface superconductivity 
have been seen in experiments on graphit o 20 ^ 21 . The en- 
hanced superconducting density has been reported on 
twin boundaries in Ba(Fei_ x Co a: )2As2 2 ^. These obser- 
vations might be explicable with our theory. Our predic- 
tions may be used for search or for artificial fabrication 
of layered and/or twinned systems with high- and even 
room-temperature superconductivity. 
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